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ABSTRACT 

Gravitational waves are messengers carrying valuable information about their sources. For sources at cos- 
mological distances, the waves will contain also the imprint left by the intervening matter. The situation is in 
close analogy with cosmic microwave photons, for which the large-scale structures the photons traverse con- 
tribute to the observed temperature anisotropies, in a process known as the integrated Sachs-Wolfe effect. We 
derive the gravitational wave counterpart of this effect for waves propagating on a Friedman-Robertson- Walker 
background with scalar perturbations. We find that the phase, frequency and amplitude of the gravitational 
waves experience Sachs-Wolfe type integrated effects, this in addition to the magnification effects on the am- 
plitude from gravitational lensing. We show that for supermassive black hole binaries, the integrated effects 
could account for measurable changes on the frequency, chirp mass and luminosity distance of the binary, thus 
unveiling the presence of inhomogeneities, and potentially dark energy, in the Universe. 



Observations of gravitational waves (GWs) have a tremen- 
dous potential for transforming our understanding of the 
cosmos. For sources at cosmological distances, such as 
the inspiral of supermassive black hole (BH) binaries in 
the sensitivity window of the proposed Laser Interferome- 
ter Space Antenna (LISA), GWs will not only carry infor- 
mation of the characteristics of the source but also contain 
information on the cosmological expansion of space-time 
through which these waves propagate, yielding redshifted 
measurements of the pro perties of the binary (e . g. distance , 
chirp-mass, frequency 
Flanagan and Hughes 
cident electromagnetic (EM) counterpart, GWs observations 
from supermassive BH inspirals might have the potential to 
serve as standard sirens for det ermining the distance-redshift 
relation dHolz and Hughesll2005l) . 

In addition to the cosmological expansion effect, GW prop- 
agation is also affected by (1) the proper motion of the source 
and the receiver relative to the cosmological flow, (2) the 
gravitational potentials at the emitting and receiving loca- 
tions, and (3) the intervening matter the GWs traverse. These 
processes are identical to those experienced by cosmic mi- 
crowave background (CMB) photons, which lead to tempera- 
ture anisotropies AT /T given by 



e pro perties of the binary (e . g. distance , 
ic v) dHughes and Hold 120031: lKnr3 [1995 
i[l998). When combined with a coin- 



= n-v\ r e + <S>\ r e + 2 f d T $dX, 



(i) 



where $ is a scalar perturbation, v is the source velocity, 
n = x/r is a unit vector, and the sub- and super-scripts r and 
e stand for evaluation at the receiver and emission location, 
respectively. In Eq. (|T), the first term is the Doppler cor- 
rectio n; the second term g ives the ordinary Sachs-Wolfe ef- 
fect dSachs and Wolfdlj.9671) : and the third is the integrated 
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Sachs - Wolfe (iSW) or Rees-Sciama effect dRees and Sciamal 
1968). The iSW term is an integral along the photon geodesic, 
with affine parameter A, which accounts for a net gain or loss 
in the temperature of the CMB photons due to changes in the 
gravitational potential of scalar inhomogeneities. This effect 
has proven to be a crucial i n interpreting observations of CMB 
inhomogeneities (see eg. dGiovann ini 2005; Afs hordil 12004) 
for a review). 

In this Letter, we present a derivation of the GW counter- 
part to the EM iS W effect. We focus on GWs propagating on a 
Friedman-Robertson- Walker (FRW) background with generic 
scalar perturbations. The d erivation follows Isaacson geomet- 
ric optics approximation dIsaacsonlll968allb]) . We find that 
the GW phase and frequency experience changes similar to 
those in Eq. ([TJ for CMB photons as the former propagates 
on a perturbed FRW background. In addition, we derive the 
corresponding changes to the GW amplitude, including mag- 
nification effects due to gravitational lensing. We estimate 
the impact of these integrated effects on measurements of the 
chirp mass, frequency and luminosity distance of supermas- 
sive BHs at cosmological distances and discuss their potential 
for unveiling the presence of inhomogeneities in the Universe. 
Latin letters from the beginning and middle of the alphabet 
will denote spacetime and spatial indices, respectively. We 
use geometric units G = c= 1 . 

Isaacs on's Geometric Op tics Approximation. Following 
Isaacson (Isaacson 1968a b), the space-time metric is decom- 
posed as g a b = "fab + chab, where 7 a j is a background metric 
and hob a GW metric perturbation, e is a book-keeping pa- 
rameter to keep track of the order of the perturbation; for- 
mally e ~ X/L <C 1, with A the wavelength of radiation and L 
the radius of curvature of the background. At the end of the 
calculation, e is set to unity. 7^ is in addition decomposed 
as "f a t = 7a& +7oj > where the labels in parenthesis denote the 
order of the perturbation. In our case, 7^ = a 2 (T)r] a b is the 

flat FRW metric and 7L = a 2 (T)5r] a b its first order perturba- 
tion. We consider only scalar perturbations (<&,$), such that 
5r)abdx a dx h = -2 [$ dr 2 + * rjijdx'dx'] . 
The linearized Einstein equations to first order in e and in 
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Lorentz-gauge (y hh ab = 0) are given by: 



V c V c h ab + 2R cadb h 



cd 



hcdT Cdr )ab + h c dl Cd ( T ab - -7 'lab 



ab 

i 



(2) 



In Eq. (0, Rgbcd, Tab and V„ denote respectively the curvature 
tensor, the stress-energy tensor and the covariant derivative 
associated with the background metric 7^. Tab is the stress- 
energy tensor associated with h a b, whose trace-reversed form 
is h a b = h a b ~ labh /2, where h and T are traces with respect to 
the background. The stress-energy tensor T a b is decomposed 
as Tab = T$> + T^ with T (0)ab = {p + p )U <0)a U (0)b + p>y ( - 0)ab 
and r (1)()0 = a- 2 (5p-2p§), r (1)() > = a\p + p)v j , and = 
a~ 2 (5p + 2pty)ri'j . Similarly, the four- velocity of the fluid in 
the background geometry is given by U a = i/ (0 - )a + t/ (1)fl with 
yds* _ (i 5 0)/a and U m " = (-$, v')/a. The quantities p and p 
are the density and pressure of the background fluid, while 8p 
and 5p are the density and pressure perturbations, where v ! is 
the 3-velocity of the perturbed fluid. 

The field equations can be simplified by choosing the TT 
gauge in this local Lorenz frame: /zoo = hoj = and hijrfi = 0. 
Thus, habl ab = and T ab h ab = 0. In addition, we neglect the 
response of the matter background to the presence of the GW 
h„b and set Tab = 0. Therefore, Eq. (0 becomes: 



V c V c hab + 2R cad bh cd = 



(3) 



Under Isaacson's (Isaacson 1968a b) shortwave or geomet- 
ric optics approximation, the GW can be written as 



hab = A a b e 



We , 



' ^ab *A ^ 



i^/e , 



: e ab h , 



(4) 



where A and <fi are real functions of retarded time u = T — r, 
e a b is a polarization tensor and r is the distance to the source. 

Isaacson's approximation allows us to simplify the field 
equations dramatically. With Eq. (0, Eq. (0 becomes 



€ 2 [-k c k c A ab ] + ie~ l [2 k c V e A fl& + A ah V c k c ] 



+[V c V c Aab + 2R cad bA cd ] = 0, 



(5) 



where k a = V a </> is the GW wave-vector. To 0(e~ 2 ), Eq. (0 
requires that k a k a = 0, implying that GW rays are null vec- 
tors and the curves x a (l), defined by dx" /dl = k a , are null 
geodesies, i.e. k b V h k a = 0. To £>(e _I ), Eq. © implies that 
k c V c e a b = 0, i.e. the polarization tensor is parallel-transported 
along null geodesies, and thus 



d_ 

dl 



\nA- 



-y a k a , 



(6) 



where we have used d/dl = k a V a . Eq. © shows that the GW 
amplitude decreases as the null rays diverge. This equations 
can also be rewritten as V ' a (A 2 k") = 0. 

Integrated Sachs-Wolfe or Rees-Sciama effect. Since the 
GW wave-vector k a satisfies the null geodesic equation, one 
can essentially follow th e derivation of th e iSW effect for 
CMB photons as given bv lPvne and Carrolll (119961) . The first 
step is to notice that the null geodesies x a {l) with affine param- 
eter I of GWs in the background metric j a b = ^(Vab+^Vab) are 
the same as the null geodesies x^A) with affine parameter A 
in the perturbed Minkowski metric 7^ = r] a b + 5riab- The affine 
parameters, metrics and wave-vectors are related by dl = adX, 
"fab = a 2 jab and k" = a k a , respectively. 



We will set coordinates such that the observer is at the end 
of the GW world-line, receiver's location jc"(r, ) = (t>,0), with 
the world-line starting at the "surface" of emission defined by 
the spacelike hypersurface of constant conformal time r e . The 
perturbed null geodesic f°(A) = i (0)fl (A)+x (1 ^(A) and its cor- 
responding wave-vector fc°(A) = k (0)a (X) + £ (1)fl (A). The lowest 
order contributions are given by x (0)fl (A) = [A,(A r - X)n'] and 
£(0)a _ (j _ n t^ ^ w here the vector n' points in the sky direction 
of arrival of the GW. 

The next order wave-vector can be obtained from the null 
geodesic equation associated with the perturbed Minkowski 
metric 7^: 



dX 

where we have used that T 



Eq. (0 yields 



(0)a 
be 



(7) 

0. The time component of 



where d$/dX = <9 T $ + fc (0) '<9, < I> and its integration yields 
fc(U0 = _($ + $)| -2$|i+/ ISW 



(8) 



(9) 



where we have defined 



i m = J d T (y+$)dX'. 



The spatial component of Eq. (0 yields 



—(*-$)+ a r ($+*) 

dX 



dX 11 

±W = -(V ij ~ k (0)r k^)dj(<P + ¥) , 
d A 



(10) 



(11) 



(12) 



where we have introduced the notation k^P' = k^'kjkP^ and 



--±\ tflV, such that k (1)i 



ki +k\ . The perpendicular 



-L ' """" " -~\\ -L 

operator _L'.= Sj—kf^kj projects components of tensors or- 
thogonal to the unperturbed wave-vector k (0)l , and thus, "par- 
allel" and "perpendicular" are operations defined with respect 
to this vector. Integration of Eqs. ( fTTT i and ( TT2T > yields 



Lav. 



(*-$)li,+/ B 



(13) 
(14) 



■± ij J dj($ + *)d\', 

where the integration constants in Eqs. (0, ( fT3l ) and ( TT4b are 
chosen so that k" is null at X e . 
The GW phase <f> is obtained from 

d 



d ^ = k { ^a^ = -k m ° + kf ) kf -. 



which after integration yields 



(* + $)£/ A'. 



(15) 



(16) 



This phase shift corresponds to the Shapiro time delay com- 
monly associated with photon propagation. 

The GW frequency lu in the reference frame of the cosmo- 
logical fluid defined by the four-velocity U a = (1 - $, v')/a is 



Gravitational Integrated Sachs-Wolfe Effect 
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given by u> = -U"k a = a" 1 [ 1 - \ e ) - $ | * + ft ■ v + I ISW ] . From 
this, one finds that the emitting and receiving frequencies are 
related via 



U e f e (1+z) 



(17) 



where \ = ^le - " ' ^"-^(-V) and a r /a e = 1 + z. Notice that 
the redshifted frequency acquires an iSW correction, as well 
as a non-integrated Doppler one. 

The GW amplitude A is obtained from Eq. ©. To lowest 
order in the scalar perturbations, we have 



4lnU ( %) = --Vf£ (0)a . 
dX 2 



(18) 



which implies that „4 (0) ar = Q is constant along the null 
geodesic with r the areal coordinate distance defined by the 
background metric rj a b. The quantity Q is determined by the 
local wave-zone source solution. At the receiving location, Q 
is given by the same solution evaluated at the retarded time. 
To next order in the scalar perturbations, Eq. (O yields 

-2 4rti = d T k mo + dM )i + dW+f a ^~k ( - 0)b , (19) 
dX 11 

where we have introduced A = ^4 (0) (1 + £) and each term is 
given by 

cU (I)0 = d T (-2$ + / ISW ) 
dik™ = - $ + 7 ISW ) - d T (V - $ + 7 ISW ) 



r™ fl p)*=— ($-3*). 

a dX 

Above, we have ignored terms of 0(r~ l ) since the l/r dom- 
inant dependence in A has been explicitly accounted at the 
lowest order. Equation ( [T9| > can be written as 

2^=-a T (*+$)+-^-(-2*+/ ISW ) 

dX dX 



-±' J I dij(<f> + y)dX', 



(20) 



which after integration yields 



$=-*l^-I±yy j d u (<!>+Wdx'. (2i) 

Notice that the iSW contribution to the amplitude has can- 
celed. The remaining integral c ontribution is the magnifica- 
tion due to gravitational lensing (Takahashi 2006). 
Combining all results, the GW takes the form 



h = Ae 



2(1 +z) 



(22) 



where we have set e to unity, 5<fi and £ are given by Eqs. ( TToT l 
and (f2TT >. respectively, and di = ar( \+z) is the luminosity dis- 
tance. At Newtonian (quadrupole) order and for an inspi- 
raling binary dPoisson and Willi [l99lh . Q = M e {-n f e M e fl 3 
and 4> e = (f} c -(-Kf e M e T 5 / 3 /16, with M. e and f e the intrin- 
sic "chirp mass" and frequency of the binary, cf> c the value of 



the phase at / = oo and t(f) = t c -(5/256)M e (irf e M e y s/3 - 
Therefore, Eq. ( 122b becomes 

h = Mei } +Z \ nf e M e ) 2 '\\ . (23) 

d L 

The modified redshift relation [Eq. ( TlTb l implies that f e M e = 

frM e (l+z)/(l-xl and thus M r = Med+z)/(l-x) and 
Eq. (l23l becomes 



(24) 



where Dl = ^z./(l — X+0 is me modified luminosity distance 
and 0,. is the modified version of Eq. ( TToT ). 

The Fourier transform of Eq. d24l >. using the station- 
ary phase approximation, is h = (JvQ/Di)(f, ■M,)~ 1 ^e'^ r , 
where we have neglected the antenna-pattern functions and 
where ip r = 2ir f r t a + <fi r (t ) with t„ a stationary point of the 
phase. The square of the signal-to-noise ratio (SNR), a 1 = 
4 J °° \h\ 2 /S„df r , is then given by 



■Dz. 



(frM r r V3 

(5„/X r ) 



dC/ r 7W r ), (25) 



with S„ the spectral noise density. 

Therefore, the changes to the chirp mass, frequency, lu- 
minosity distance and SNR induced by the scalar inhomo- 
geneities in the background are given by 



SM _ Sf _ 
M z ~ f z ~ X ' 



-r- = x-£, — = t + ^' (26) 

di a L 



where M z = M e (l +z) and f z = f e /(l +z). 

Root-Mean-Squared Fluctuations. Next, we expand \ and 
£ in Fourier modes to relate them with the density perturba- 
tion modes: = = -3H$Sf' 'D(t) /(2k 2 a), where D(?) is the 

linear growth function and <5^ 0) is the Fourier coefficient of the 
density perturbation at zero redshift. The ensemble average is 
performed over the density perturbations via the definition of 
the power spectrum 



S^St^) = 5\k-k')(2n) 3 Pf\k) = 



25^ 



PrT z 



(27) 



where P R (k) = 2tt 2 A^ 0) V 3 is the power spectrum of curvature 

perturbations, with A® = 2.445 x 10" 9 , and T(k) is a transfer 
function, which we approximate via the fitting relations given 
in (Eis enstein and Hulll997l) . The e' k ' x factor that arises in the 
Fourier transform is expanded in spherical Bessel functions 
and the integration over th ese functions is performed th rough 
the Limber approximation dLo Verde and Af shordi 20081) 



dt'—^jJ i+l / 2 {kt') • 



x/W+T/2) 



f 



1+1/2 



(28) 



which requires that t+ 1/2 = kt ', valid for large I. 
With these considerations in mind, 




kdkP R 



k 5 dkP R 



1/2 



1/2 



dt'g\t') 



1/2 



(29) 



dt't' 2 g 2 (t') 



1/2 

(30) 
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FIG. 1. — Root-mean-squared fluctuation of \ (dotted) and £ (dashed) as a 
function of redshift. 



where k min « 0.0004, g(t) = D(t)/a(t) and Cl mfl « 0.27. 
The g is approximated via the fitting function 
D/D = nl ^q-^H- 2 ^ 1 (lAcauaviva et all l200l 
iPolarski and Gannouiil [2008b . where j « 0.557 - 0.002z, 
H is the Hubble parameter, and Hq is its value today. a(t) 
and H(t) are obtained by solving the Friedman equations 
nu merically for a cos mology with parameters measured 
in dKomatsu et al.ll2009ft . 

FigureQ]shows the root-mean-squared fluctuation of x (dot- 
ted) and £ (dashed) as a function of redshift. The vertical line 
approximately corresponds to the limit to which the Advanced 
Laser Interferometer Gravitational Observatory (LIGO) will 
be able to see for an optimally-oriented binary. We have here 

1/2 

somewhat underestimated the (x 1 ) term, as have neglected 
non-linear corrections to the power spectrum d ue to galaxy 
and cluster formation (Peacock and Dodds 1996). 



Data Analysis Implications. The detection of GWs by LISA 
or any other instrument will not be affected by matter in- 
homogeneities, since the phase correction accounts for a to- 
tal phase shift that is extremized over during GW extraction. 
Parameter estimation, however, will clearly be affected by 
matter inhomogeneities. LISA is expected to be sensitive 
to the chirp mass and the luminosity distance to Aln.M ~ 
10" 6 and AlnD L - 10" 3 dArun et al J 120071) for low-redshift 
sources (z ~ 0.55), but also sensitive to high redshift sources 
(z ~ 5 - 10) up to AlnD L - 10" 2 dHughesI 12001 IVecchiol 
2004). Our calculations confirm that £ (associated with weak- 
lensing) will be a noise source in the use of standard sirens 
to measure the equation of sta te of dark energy through the 
redshifted luminosity distance dKocsis et al.ll2006l) . 

Alternatively, one could view these effects as a new link 
between EM measurements of density inhomogeneities and 
LISA observations. In order to achieve this goal, however, one 
would first have to break the degeneracy between (M. z ,z,X) 
or (cIl,z,£,,x)- Given a coincident EM and GW detection, 
one might be able to achieve just this, by electromagnetically 
determining the redshift and the component masses via host 
galaxy identification and correlations between galaxy lumi- 
nosity and BH mass. Another possibility would be to use 
large scale structure observations to measure 8k and predict x- 
Cosmologists already use large scale structure obser vations 
to predict the iSW term for the CMB (|Afshordili2004|)- Such 
measurements would then open up, for the first time, studies 
of cross-correlations bet ween GWs and large-scale structure 
surveys of dark matter (Crittenden and Turok 1996) and pos- 
sibly dark energy dCooray et al.ll2004l) . Furthermore, a detec- 
tion of the cross-correlation between matter distribution and 
the GW iSW effect could potentially be another test of GR 
since it would show that GWs propagate in the same metric 
as EM radiation. 
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